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Abstract 

For a quantum mechanical system with broken super symmetry, we present 
a simple method of determining the ground state when the corresponding 
energy eigenvalue is sufficiently small. A concise formula is derived for the 
approximate ground state energy in an associated, well-separated, asymmetric 
double-well-type potential. Our discussion is also relevant for the analysis of 
the fermion bound state in the kink-antikink scalar background. 
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Supersymmetry(SUSY) and its breaking are fundamental issues in theoretical particle 
physics. There have also been numerous applications of SUSY to quantum-mechanical 
potential problems [1,2], based on the observation that the spectrum of the Hamiltonian 

H + = ~ + V + (x), V + (x) = W 2 (x) + W'(x) (1) 

(W(x) is the superpotential, and we set h = 2m = 1) is related through SUSY to that of 
the partner Hamiltonian 

d 2 

H = -— + V-(x) , V-(x) = W 2 (x) - W'(x) . (2) 

This formalism has provided us with a number of exactly solvable quantum mechanical 
systems for which energy eigenvalues and eigenfunctions can be found in closed forms. The 
key properties that made such feat possible are unbroken SUSY, manifested by the vanishing 
energy for the ground state of H- (or H + ), and shape invariance under the change of 
parameters for the given potentials [3]- [5]. This approach can sometimes be extended 
to parameter ranges corresponding to (spontaneously) broken SUSY [6,7]. But, with SUSY 
broken, the ground state energy is no longer equal to zero and this jeopardizes the possibility 
of obtaining exact analytic results by the SUSY-based method in a crucial way. 

In this work we will show that, in some broken SUSY case for which the lowest energy 
E(> 0) for the Hamiltonian H + or H is sufficiently small, a simple perturbative scheme 
leading to an easy evaluation of E can be developed. Our discussion finds useful application 
in studying the almost-zero-energy fermion modes in the background of a soliton-antisoliton 
pair. 

The superpotential relevant for our discussion is given as follows. Let ctr(x) be a generic 
function with the properties 

(7r(x) > , for x > and \x\ not very small , 
o~r(x) — > —v , for x < and \x\ not very small 

and (Tl{x) the one with the properties 

a L (x) > , for x < and \x\ not very small , 
—> —v j for x > and \x\ not very small 

so that the related potentials Vr±(x) = a\ ± o' R {x) and Vl±(x) = o\ ± o' L {x) may have the 
general shapes shown in Figs.l and 2, respectively. Then the superpotential appropriate to 
our case is obtained by combining these two types of functions as 

W(x) = a R (x - h) + a L (x - l 2 ) + v (5) 

with L = \li — l 2 \ taken to be reasonably large (so that W(x) may have a flat basin between 
the points x — l 2 and x = h). 1 See Fig. 3 for the schematic plots of W(x) and the related 



l Kt the later stage of our discussion we will use the fact if /* denotes a certain point in the fiat 
basin, the approximation <jr(x — /i) + v ~ for x < I* or a^ix — h) + v ~ for x > I* is valid. 



2 



potentials V±(x) = W 2 (x) ± W'(x). Both W(oo) and W(— oo) being positive, this corre- 
sponds to the case of broken SUSY [1,2]; but, for the posited superpotential (with L large), 
the ground state energy E is expected to be rather small. [Our superpotential will be an 
even function of x if (Tl(x) happens to be the mirror image of ctr(x), i.e., (Jl{x) = <tr{—x), 
and (h, l 2 ) are equal to (|f, — §)]. 

For W(x) specified as above, the corresponding Hamiltonians H± involve the potentials 
which can approximately be described by the sum of two well-separated potentials (aside 
from a constant term), i.e., 



V±{x) ~ V R± (x - h) + V L± (x - l 2 ) - v 2 



(6) 



These correspond to asymmetric double wells even when W(x) is an even function, and hence 
the well-known approximation schemes used for symmetric double wells (e.g., instanton 
methods, tight-binding approximations) would not be much useful. [Note that, for the 
ground state of our Hamiltonian with the potential V + ~ Vr + +Vl + — v 2 (or V- ~ Vr-+Vl_ — 
v 2 ), the tight-binding approximation is plainly not available — while the local Hamiltonian 
involving Vl+ (or Vr_) allows a zero-energy bound state, no zero-energy state exists for 
the local Hamiltonian involving Vr + (or Vl_)]. But, for a supersymmetric system, one can 
always consider a pair of coupled first-order differential equations instead of the second- 
order Schrodinger equations. Our perturbative approach for the ground state is based on 
the analysis of these first-order equations, and as a result we obtain a remarkably simple 
formula for the lowest eigenvalue E. It is simply the square of the product of the two zero- 
energy eigenfunctions (allowed with the potentials Vl + (x — l 2 ) and Vr_(x — l\) separately), 
evaluated at an arbitrary chosen point I* in the flat middle region of the superpotential 
W(x). See the expression (21) below. 
Consider a matrix Hamiltonian 

with 

A = d x + W(x) , A ] = -d x + W{x) . (8) 

f^i(x)\ 

The corresponding eigenvalue equation, 7i^{x) = u^/{x) with ty(x) = T ; (), reduces to 



a pair of first-order differential equations 



*2(X) 



AV 1 (x)=uV 2 (x) , (9a) 

A ] V 2 (x) =wVi(x) . (9b) 

Then, based on the relations 

A^A = —8% + W 2 (x) - W'(x) = H , AA^ = -d 2 x + W 2 (x) + W\x) = H + , (10) 

one finds that the functions ^i(x) and ^ 2 (x) are eigenfunctions of the Schrodinger Hamilto- 
nians H_ and H + (with the same energy E = u 2 ), respectively. Further, it is not difficult to 

f^i(x)\ ~ ( ^i(x) \ 

show that if ^(x) — T ; ( is an eigenvector of H with eigenvalue u, ^f(x) = T , \] 
v ; U 2 x i V ; \-y 2 (x)J 
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corresponds to an eigenvector of H with eigenvalue — oo. For our later application, we also 
remark that the above matrix Hamiltonian can be written as 



U = -i^^d x + -y°W(x) 



(11) 



with 7 = a 1 = 



1 

1 



and 7 1 = icr 3 = 



i 
-i 



. That is, it describes the Dirac Hamiltonian 



(defined on a line) in the presence of an external scalar field W(x) [8]. 

Let us now specialize to the case with the superpotential given by the form (5) with L 
large. If tpi(x) is the ground-state eigenfunction of the Hamiltonian H_ with a small energy 
E = u) 2 and <f 2 (x) that of the isospectral partner H + with the same energy, they will have 
to satisfy the integral equations (derived from Eqs. (9a) and (9b)) 



Pi 



(x) = e~^ W{y)dy {D l +u 



dy 



f, v W(z)dz / x 



ip 2 {x) 



e Ii X 2 W(y)dy 



Do — 



CO 



dy 



OL2 



■fi y W{z)dz 



<pi{y) 



(12a) 
(12b) 



for suitable constants a±, a 2 , Di and D 2 . Clearly, for small enough u>, we may iterate these 
equations to learn about the corresponding eigenfunctions. Then, let (Pi(x) and (p 2 (x) denote 
the limiting expressions of <fii(x) and (f 2 (x) for very large L (and therefore very small u). 
Here, considering the general shapes of the potentials V± (see also Eq.(6)) shown in Fig.3, 
(fi(x) {p 2 {x)) can be taken as the normalized zero-energy eigenfunction of — d 2 + Vn-(x — 1\) 
{—d 2 + Vl + {x — l 2 )). This in turn implies that <fi(x) and (f 2 (x) satisfy the first-order 
differential equations 



[d x + <tr(x - h)]ip\{x) = 
Hence we have the explicit expressions 

( p° 1 (x) = C 1 e- I "i <TR{v - ll)dv 



[-d x + a L (x - l 2 )]tp° 2 (x) = . 



tplix) 



C 2 e Il 2 aL{y ~ h)dy 



(13) 



(14) 



with appropriate normalization constants C\ and C 2 , which are taken to be positive. [In view 
of Eqs. (9a) and (9b), both (fi(x) and <p 2 (x) are chosen to be real]. As these informations are 
used with our integral equations (12a) and (12b), we are led to the following approximate 
expressions for <fii(x) and tp 2 (x) (valid for large L): 



ipi{x) 



p 2 (x) 



) ~ e~ J"'" W(y)dy \Ci+Q dy W{z)dz ?l{y) 



^ e Ii x 2 w(y)dy 



Co. 



dy 



■ III W(z)dz r n 



(15a) 
(15b) 



with the value I* chosen conveniently at some point in the flat middle region of the super- 
potential so that both | Zx — l*\ and \l* — l 2 \ may become 0{L). 

We then note that the expression (15a) as x — > — oo (and similarly that in Eq.(15b) as 
x — > oo) would blow up unless the value of uj were chosen such that 



Ci 



dy 



W(z)dz 0/ \ 



. 



(16) 
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Here, with the explicit expressions for W(x) and ^"(x) given in Eqs. (5) and (14), we observe 
that 



e J l y l W(z)dz _ ^ a R {z-h)dz e ^a L {z-l 2 )dz e !^ L (z-h)dz 

= CMin^r)}- 1 ^) , (for ye [-oo, /•]) 



(17) 



Ji* 



/oo 
dvWv)? = -1 , (18) 
-oo 

and hence the 'eigenvalue condition' (16) reduces to the form 

c, - ioc^in^n]- 1 = o . (19) 

We thus find 

u = <p1(l*)<p%(l*) = dC 2 e~ S h '"(v-Wy+fh °l(v-Wv (20) 
and so, for the desired ground-state energy, the formula 

e = [<A(n<fi(n] 2 • (2i) 

[Notice that, for an arbitrary choice of I* in the flat middle region of the superpotential, 
the same value for E results]. The same expression for E is obtained by analogous analyses 
with Eq.(15b). With the eigenvalue u determined in this manner, the corresponding wave 
functions ipi(x) and ip?.{x) (up to normalization) are now expressed as 

(u[^(x)]-^-dy[^(x)]\ x>l* 

^ ] ^^){i+ w /;«rM(!/)}, x<l*. { } 

We expect that a judicious use of the tight-binding approximation with the matrix Hamil- 

fm ( X )\ ( \ 

tonian (7), taking ( q J an d ( (p°( x ) ) aS ^ e ^ e S enera ^ e (*' e -> zero ener gy) eigenstates of 

the corresponding local Hamiltonians, lead to the same conclusion as above. 2 This is sup- 
ported by the observation that, for the eigenvalue u>, the same result (i.e., Eq.(20)) follows 
from the calculation based on the formula 

/oo 
dx <pl(x)H(p±(x) , (24) 
■oo 

/ m°(x) \ 

taking tp±(x) = ( _|^o^ )> the usual zeroth-order states in the tight-binding approxi- 
mation. Indeed, one can verify that the given integral 

/oo poo 
dx ip + (x)H<p° + (x) = dx [a R (x - h) + v\y\{x)yl{x) (25) 
-oo J — oo 



2 But, since our matrix Hamiltonian (which is the Dirac Hamiltonian) is unbounded from below, 
some care must be exercised. 
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is well approximated by the product tp^l*)^^*)- 

As an explicit example, consider the superpotential of the form 

W(x) = v sgn(x — — ) — v sgn(x + — ) + v , 



(26) 



i.e., in our notation, (I1J2) = (§,—§), cr R (x — h) — v sgn(x — ^), and a L (x — l 2 ) 
—v sgn(x + §)• Given this, the potentials of the Schrodinger Hamiltonians H± will be 



V±{x) =v 2 ±2v 5(x ) =f 2v 5(x + -) . 



(27) 



For these systems, one can of course find the exact ground state energy E by solving the 
appropriate Schrodinger equations. This exercise shows that E is the root of the equation 
E = v 2 e - 2L V v2 -E ; and hence, for large L, we have 



E 



~ v 2 e~ 2vL 



(28) 



Let us see whether our formula (21) yields the same. The normalized solutions of Eq.(13) 
are trivially found here: 



v e 



M*-k\ 



<p° 2 (x) 



v e 



-v\ x +k\ 



(29) 



Then, from Eq.(21), we have that E (for large L) should equal [</??(0)<^(0)] 2 = v 2 e~ 2vL . 
Hence a complete agreement. 

More physically relevant example is provided by the Dirac Hamiltonian (11) with the 
scalar field taken to represent the kink-antikink pair, 



W(x) = vtanh 



v tanh 



A 4 / L 
-(x + - 
2 V T 



+ v . 



(30) 



Here the scalar field <Jr(x — ■§ 
and <Tl(x + ^) = 



'tanh 



^)] represents a kink located at x = 
tanh [^(x + ^)] an antikink at x = — |f [9]. The widely-separated 
kink-antikink configuration, described by the form (30), has received attention in Refs. 
[10-12]. Especially interesting is the almost zero-energy mode of the Dirac Hamiltonian, in 
connection with the role of the so-called Jackiw-Rebbi mode [10] (which refers to the zero- 
energy fermion mode [13] in the kink or antikink background) when a kink and an antikink 
are simultaneously present. In the kink or antikink background the Jackiw-Rebbi mode is 
V?(x)\ ( 



represented by 
Eq.(13): 







or 



, if <Pi(x) and tp\ (x) denote the normalized solutions of 



cosh I (x- ^) 



/in 



<P°{x) = C 



, A 4 / L. 
cosh-(x + -) 



(31) 



where C 



p. J* 2 ' 
2r(f )r(|) 



. Then, in the above kink-antikink scalar background, one can 



immediately find the energy of the almost-zero-energy fermion eigenmode by using our 
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formula (20) — it is equal to ±o>, with 3 



u, = ^?(0)^(0) = C 2 [cosh^]-f 

rss —p u — p 

"2 e r(f)r(i) e • 

The exponential dependence of uj on the distance L was previously noted in Ref. [10]. This 
Independent fermion energy shift will contribute to the effective potential between the kink 
and the antikink. For instance, in the vacuum sector where all negative energy fermion 
modes are to be occupied, the contribution from this mode, i.e., that with energy —uj will 
become more negative as L decreases, thus producing an attractive interaction between the 
kink and the antikink (if only this mode is taken into account). 

In this work we investigated on some special properties pertaining to the ground state 
of a quantum mechanical Hamiltonian with broken supersymmetry, when the corresponding 
eigenvalue is small. A direct perturbative analysis based on the first-order coupled equaitons 
has been used to obtain a very simple (approximate) expression for the ground state energy. 
Our formula (21) should be useful in finding the ground state energy of a Schrodinger 
Hamiltonian the potential in which can be approximated by the form (6). 
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3 Based on some heuristic arguments the author of Ref. [12] also identified uj with the expectation 
value in Eq.(24). But his final expression for uj is apparently not consistent with our result in 
Eq.(32) — this is due to the calculational error made in Ref. [12] (e.g., wrong fermion zero-energy 
wave functions used). 
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FIGURES 




FIG. 2. Schematic plots of (Tl(x) and Vl±(x) = cr 2 L (x) ± <j' l {x) 



